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These entire analytic functions enjoy properties, such as rearrangement and unique 
determination by infinite jets, that one would expect. However, not all exponential 
functions (eigenfunctions of the Laplacian) arc entire analytic, and also many other 
natural candidates, such as the heat kernel, do not belong to this class. Neverthe- 
less, we are able to use spectral decimation to study exponentials, and in particular 
to create exponentially decaying functions for negative eigenvalues. 

§1. Introduction. 

Ordinary calculus is such a remarkable subject because it combines both a gen- 
eral conceptual framework and a detailed understanding of basic functions. For 
example, the theory of power series expansions hinges on the elementary observa- 
tion that the function fnix) = x^/nl on [0, 1] is bounded by 1/nl. (Stated this way, 
it seems almost a tautology, so perhaps it is better to say that fn is the polyno- 
mial characterized by the conditions fn^\o) = Snm-) Another example: among all 
linear combinations of coshx and sinha; there is one, = coshx — sinhx, that 
decays as a; — > oo; moreover its rate of decay is the reciprocal of the growth rate of 
cosh a: and sinha;. 

The goal of this paper is to understand analogous facts about basic functions 
on the Sierpinski gasket (SG), which should be regarded as the simplest nontrivial 
example of a fractal supporting a theory of differential calculus based on a Lapla- 
cian. Standard references are the books of Barlow [Ba] and Kigami [Ki2] , and the 
expository paper [S2]. The references to this paper, and the more extensive bibli- 
ography in [Ki2], indicate an intensive development of the subject since Kigami's 
original paper [Kil] giving a direct analytic definition of the Laplacian on SG. 

Recall that SG is the attractor of the iterated functions system (IFS) consisting 
of three contractions in the plane Fi{x) = ^{x + qi), i = 0, 1,2 where qi are the 

2 

vertices of an equilateral triangle. In other words = Fi{SG), and we refer to 

i=0 

the sets Fi{SG) as cells of order 1. More generally, we write = F^^ o ■ ■ ■ o F^^^ 
for a word w = (wi, . . . , Wm) of length \w\ = m, each Wj = 0, 1 or 2, and call 
Fiu{SG) a cell of level m. We regard SG as the limit of a sequence of graphs 
(with vertices Vm and edge relation x ~m y) defined inductively as follows: Fq is 

2 

the complete graph on Vb = {qo, 9i, 92}, and Vm = FiVm-i with x ~rn y i^x and 

i=0 

y belong to the same cell of level m. Then K = |^ Vmi the set of all vertices, the 

m=l 

analog of the dyadic points in the unit interval, is dense in SG. We consider Vq the 
set of boundary points of SG, and \ Vq is the set of junction points. Note that 
every junction point in Vm has exactly 4 neighbors in the graph F^- The graph 
Laplacian on F^ is defined by 

(1.1) AmU{x) = J2 (^(2/) - «(^)) for xeVm\ Vq. 

y^mX 
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The Laplacian A on SG is defined as the renormahzed hmit 



(1.2) 



Au{x) = Urn -5'^Amu{x). 



m— >oo 2 



More precisely, u e domA and Au = f means u and / are continuous functions 
and the hmit on the right side of (1.2) converges to / uniformly on Kk \ Vq. The 
Laplacian plays the role of the second derivative on the unit interval (although it 
is shown in [S] that it does not behave like a second order operator) . Thus we will 
define a polynomial P to be any solution of A-'P = for some j. More precisely, 
if we let Tij denote the space of solutions of A^~^^u = 0, then TCj is a space of 
dimension 3j + 3, and it has an "easy" basis {fnk} for < n < j and k = 0, 1, 2 
characterized by 

(1-3) A^fnkiqk') = ^in^kk'- 

In [SU] a different basis was constructed in order to develop a theory of splines. 
Here we will consider yet another basis, implicitly used in [S3] in conjunction with 
Taylor expansions, to define power series. 

The Laplacian is basically an interior operator, as (1.2) is not defined at the 
boundary (although Au = f makes sense at boundary points by continuity). There 
are also boundary derivatives. The normal derivative 



(cyclic notation qj^^ = qj) exists for every u e domA and plays a crucial role in 
the theory, especially in the analog of the Gauss-Green theorem: 



Here jj, is the natural probability measure that assigns weight 3~"^ to each cell of 
order m. The normal derivative may be localized to boundary points of any cell, 
and there is also a localized version of (1.5). At a junction point there are two 
different normal derivatives with respect to the cells on either side. For u G dom A 
we have the matching condition that the two normal derivatives sum to zero. This 
leads to the gluing property: if u and / are continuous functions and Au — f on 
each cell of order m (meaning A{uo F^) = 5~"^/ o for all words w of length m), 
then Au = / on SG if and only if the matching conditions hold at every junction 
point in Vm- 

There are also tangential derivatives 



(1.4) 




(1.5) 




(1.6) 



dTu{qj)= hm 5^{u{F^q^+,)-u{F^qj-^)) 
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that exist if « e dom A, and may be localized to boundary points of cells. In this 
case there are no matching conditions for u e dom A. However, we will show in 
Section 5 that there are matching conditions involving infinite series of tangential 
and normal derivatives valid for polynomials and analytic functions. Tangential 
derivatives were introduced in [S3]. Their true sigificance is still somewhat elusive. 
In this paper we will show that for polynomials and analytic functions the sum of 
the tangential derivatives over the three boundary points of any cell must vanish. 
In [S3] and [T2] the idea of creating a gradient of a function out of the normal and 
tangential derivatives is discussed. Here we will extend this to the idea of a jet. For 
simplicity we deal with a boundary point g^,but the definition can be localized to 
boundary points of any cell. 

Definition 1.1 : For u G dom A"^, the n-jet of u at qi is the (3n + 3)-tuple of values 
{A^u{qi).i dnA^u{q(), dTA-'u{qe)) for < j < n. For u E dom A°°, the jet of u at qi 
is the infinite set of the same values for all j >0. 

(£) 

Fix a boundary point q£. We define polynomials P^fJ by requiring that the 
j-jet at qi vanish except for one term, A^Pjf'{q£) = 1, dnA^Pj2{Q(.) = 1 and 
drPj^iQi) = I5 respectively. We refer to these functions as monomials. It is clear 

that the monomials Pj^J for < j < n form a basis of Hn- It is shown in [S] 
that they exhibit a prescribed decay rate in neighborhoods of q^, but the estimates 
established there were not uniform in j. The first goal of this paper is to obtain 
sharp estimates for ||-Pjf'* ||oo- For pjf-* and Pj^^ we prove decay estimates faster 
than any exponential. For Pj2 the situation is diff'erent; we prove an exponential 
decay of order A^"' for the specific value A2 equal to the second nonzero Neumann 

eigenvalue. This result is sharp. In fact we show that {—X2)-'Pj2 converges to a 
certain A2-eigenfunction of A. This result has no analog in ordinary calculus. 

We define a power series about q^ as an infinite linear combination of the mono- 
mials Pj^ with coefficients {cjk}- We find growth conditions on the coefficients to 
guarantee convergence. We study the rearrangement problem: given a convergent 
power seres about one boundary point, does the function also have a convergent 
power series about the other boundary points? Surprisingly, we find that it is nec- 
essary to assume a stronger growth restriction on the coefficients in order for this 
to be the case, namely 

(1.7) \cjk\ = 0{W) for some R < A2. 

We end up defining an entire analytic function to be a function represented by 
a power series with coefficients satisfying (1.7). We then prove rearrangement is 
possible at all boundary points, and in fact local power series expansions exist on 
all cells, with the estimate (1.7) preserved (in fact the same R value). This choice 
of definition means that there are some convergent power series that do not yield 
analytic functions. It also means that eigenfunctions of the Laplacian cannot be 
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entire analytic functions unless the eigenvalue satisfies |A| < A2. On the other hand 
it is easy to see that there are A2-eigenfunctions that cannot be represented by 
convergent power series, so the definition seems to be close to best possible. We 
then are able to characterize the class of entire analytic functions in domA°° by 
the growth conditions 

(1.8) l|A^«l|oo = 0{R^) for some i? < A2 

(one could also use norms). 

Our definition of entire analytic function means that a basic principle of unique 
analytic continuation holds. If we have a function defined on a cell and satisfying 

(1.8) there, it has a unique extension to an entire analytic function on the whole 
space. In fact its jet at any boundary point of the cell satisfies (1.7), and uniquely 
determines the function. This implies that a nonzero entire analytic function cannot 
vanish to infinite order at any junction point. We could also define local analytic 
functions on a cell of order m by relaxing the condition i? < A2 in (1.7) and (1.8) 
for R < 5'"A2. One could hope to have a notion of analytic continuation that would 
allow such local analytic functions to extend to larger domains. However, we have 
not been able to find any interesting examples, so we will not pursue the matter 
here. 

It is easy to extend the notion of entire analytic function to infinite blow-ups of 
SG ([SI], [Tl]). The simplest of these is 

00 

(1.9) SGo, = U F^^iSG), 

but more generally we could consider 

00 

(1-10) \]Pi"Pn---FL\SG) 

n=l 

for any choice of ji, j2, js, ■ ■ ■ ■ A function on SG satisfying (1.8) for all > 
extends to an entire analytic function on any blow-up (1.10). It is not clear at 
present which, if any, of these functions will come to play the role of special functions 
(hypergeometric, Bessel functions, etc.) in real analysis. On the other hand it is 
very easy to construct many such functions simply by taking a power series with 
bounded or sub-exponential growing coefficients. The negative results of [BST] 
mean that none of these spaces of analytic functions is closed under multiplication, 
so this precludes using many standard techniques for ordinary power series. 

Although none of the eigenfunctions of the Laplacian are entire analytic functions 
on the blow-ups, it is still important to understand their global behavior. In Section 
6 we study this problem for the simplest example SG^o and negative eigenvalues. 
It is easy enough to define the analogs of the functions cosh \/\x and sinh \/\x. In 
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fact there are three, which we call Cx{x), Sx{x) and Qx{x), characterized among 
(— A) -eigenf unctions by their 0-jet at qo, or equivalently by power series involving 
just Pji\ Pj2i ^jz' terms, respectively. The power series for C\{x) and Q\{x) 
converge on all of SG oo, while the power series for S\{x) is only convergent on 
a neighborhood of go (depending on A). Fortunately, there is another method 
available to study these eigenfunctions, called spectral decimation ([FS], [DSV], 
[Tl]). Using this method we are able to show that they exhibit an exponential 
growth as a; — > oo (or as A — > oo), and there is one linear combination, E\{x) = 
Cx{x) — Sx{x) for the appropriate normalization, that decays as x — >^ oo at the 
reciprocal rate. Thus Ex{x) is the analog of e~^^. It is not clear if there is any 
analog of e^^. 

Although we do not use power series in our study of properties of eigenfunctions, 
we can turn the tables and use facts about eigenfunctions to obtain information 
about power series. In particular, we are able to construct specific power series that 
are divergent, or power series that are convergent but not rearrangeable. We can 
also give an explanation for why the recursion relations for the size of monomials 
are unstable. 

It is interesting to speculate on possible future extensions and developments 
of our results. It is important to understand all eigenfunctions, including those 
with positive eigenvalues, on all blow-ups (1.10). There should be some sort of 
Liouville-type theorem precluding nonconstant bounded entire analytic functions 
on blow-ups without boundary. 

What is the behavior of an entire analytic function in a neighborhood of a generic 
point? Is there any notion of power series there? Are there interesting examples 
of local analytic functions with a natural domain that is not just a single cell? Is 
there a meaningful notion of analytic functions on fractafolds based on SG [S4]? 

We have seen that there is no restriction on the jet of an analytic function other 
than the growth condition (1.7). For the larger class dom A°°, is there an analog of 
Borel's theorem that an arbitrary jet may be specified at one (or all three) boundary 
points? 

In [OSY], the structure of level sets of harmonic functions on SG was elucidated, 
with the remark that certain eigenfunctions of the Laplacian have level sets of an 
entirely different nature. It is clear now that these eigenfunctions are not analytic, 
so it is reasonable to ask if anything interesting can be said about level sets of entire 
analytic functions. Another remark from that paper is that harmonic functions 
enjoy a principle called "geography is destiny." Roughly speaking, this says that 
the restriction to a small cell of a harmonic function is essentially dictated (up 
to two parameters) by the location of the cell, rather than the specific harmonic 
function, in a certain generic sense. This holds because restrictions of harmonic 
functions are governed by long products of matrices, so the theory of products 
of random matrices makes generic predictions. For analytic functions, there is a 
similar description of the transformation of jets, except that the matrices are now 
infinite. So if we go to a small cell, while all jets satisfying (1.7) are possible, some 
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may be very unlikely for a generic analytic function. Is there some way to make 
this precise? 

A sequel to this paper, [BSSY], will discuss functions with point singularities, 
exponential functions on general blow-ups, and estimates for normal derivatives of 
Dirichlet eigenfunctions and heat kernels. 
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§2. Polynomials. 

The space Tij of (j + l)-harmonic functions (solutions of A-'+^w = 0) has di- 
mension 3(^ + 1) and plays the role of the space of polynomials of degree at most 
2j + 1 on the unit interval. Several different bases for T-Lj are known. In [SU] , in 
order to develop a theory of spline spaces, bases based on the behavior at all three 
boundary points were used. In this section we will discuss properties of yet another 
basis, based on the behavior at a single boundary point, that is more suited to the 
work on power series to follow. The polynomials in this basis are analogous to the 
monomials /n\ on the unit interval. These functions were introduced in [S3] , but 
not much was done there to describe their behavior. 

Definition 2.1 : Fix a boundary point q^. The monomials PS for k = 1,2,3 and 



J = 0, 1, 2, . . . are defined to be the functions in Tij satisfying 

(2.1) A^P^'^{qi) = SmjSki 

(2.2) dnA'^P^'^{qi)=6mjSk2 

(2.3) dTA^P^l\qe)^SmjSks- 



When £ = we will sometimes delete the upper exponent and just write Pjk- 

Note that we only need to consider m < j in (2.1-3), since A"^Pj^ vanishes 
identically otherwise. Thus there are 3(j ' + 1) conditions in all, and it follows from 
[S3] that there is a unique solution, and the monomials Pj^ for fixed £ and all 
j < ji form a basis for Tij^^ . We have the self-similar identities 

(2.4) P^f\Fpx) = 5-^-^P^i\x) 

(2.5) P^2\Frx) = (ly^-^-^P^^ix) 

(2.6) P^SiFtx) = 5-^^+'^"^P^i\x) 

that describe the decay rate of these functions as x ^ q£ (of course Pq^ = !)• It is 
easy to see that P^f^ and P^^ are symmetric while Pj^^ is skew-symmetric under 
the reflection that flxes qi and permutes the other two boundary points. It is easy 
to compute the values of monomials to any desired precision. Figure 2.1 shows the 
graphs of some of them. Since we may obtain P^^ from Pj'^^ by simply rotating 
the variable x, we will restrict our discussion to £ = from now on. 
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Figure 2.1 : The graphs of Pjk for some typical values. The graphs of Pji are all 
qualitatively similar for j > 1, so we show only P51. Similarly for Pj^- The nature 
of the graphs of Pj2 changes drastically around j — 5, 6, 7, 8, so we display all of 
these. The graphs of Pj2 for j > 8 are qualitatively similar to P82- 
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It is clear from the definition that powers of the Laplacian send monomials to 
monomials, simply reducing the j index: 

(2.7) ^"^Pjk — P{j-m)k- 

We could use this property to give an inductive definition. When j = the 
monomials are explicit harmonic functions, Pqi = 1; Pq2 has boundary values 
^^02(<?o) = 0, Po2(9i) = ^^02(92) = -1/2 and Pqs has boundary values Po3(?o) = 0, 
-Fo3('?i) = —Posi(l2) = 1/2. Then Pj^ for j > is the unique solution of APj/j = 
P{j-i)k with vanishing initial conditions 

Pjk{qo) = 0, dnPjkiqo) = 0, drPjkiqo) = 0. 

In [KSS] it is shown that Pjk may then be written as an integral operator (with 
explicit kernel) applied to P(j_i)/i;. However, the kernel is quite singular, so we have 
not been able to extract any useful information out of this representation. 

There are three main goals in this section: 1) to obtain sharp estimates for the 
size of the monomials, 2) to understand how to express monomials for one choice 
of £ in terms of monomials for another choice of 3) to obtain certain universal 
identities that hold for all monomials. In pursuit of these goals we introduce some 
terminology. 

Definition 2.2 : For j >0 let 

\ Hj = dnPjl{qi), tj = dTPj2{qi)- 

Note that by symmetry we have Pji{q2) = otj, Pj2{q2) = I3j and Pj^{q2) = — 7j, 
so that all values of monomials at boundary points are expressible in terms of a's, 
/3's and 7's. Soon we will see that the n's, t's and ct's suffice to express all normal 
and tangential derivatives of monomials at boundary points. 

Theorem 2.3. The following recursion relations hold: 

4 

(2.9) aj = — — -J^oij-iOie for J > 2 

^ ^ ^ 1=1 



(2.10) 



7i 



5^+1 - 5 



(2.11) p. = -±-J2[l5^-%-tPi-laj-i5'Pt + laj-iPt) for j > 1, 
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with initial data ckq = 1, cti = 1/6, /3o = — 1/2, 70 = 1/2. In particular, 
(2.12) 7j = 3aj+i. 

Proof : It is convenient to work in matrix notation, with all matrices being infinite 
semi-circulant. For example, the matrix a = {aij}ij=o^i^2,... has n for 

i > j and = for i < j. We consider two linear operators on such matrices, the 
shift a and the dilation r, given by 



T 



/do ■■■ 

di do 

d2 di do 



fdo ••• 

di do 

d2 di do 



J 



^rfi ■■■ ^^ 

d2 di 

ds d2 di 



V 



\ /do ••• \ 

5di do 
52^2 5di do 



Let {/ji, fj2, fj3}'jLo be the easy basis defined by (1.3). As in [SU] we let 

ai-i = dnfikiqk) 

bi-i = dnfikiqn) n^k 

for Z = 0, 1, 2, . . . . Then the Gauss-Green formula says for Z > 

m = dnf(i+i)i{qi) 

3 

= Yl (/oi(5n)5n/(/+l)l(gn) - (gn)9n/oi (Qti)) 

n=l 

= / (/oiA/(/+i)i - /(/+i)iA/oi)cZ// 

JSG 



and 



= / foifiid/J, 

JSG 

k = dnf{i+i)i{q2) 
3 

= if02{qn)dnf(l+l)l{qn) - /(Z+l)l (?n)5n/o2(gn)) 

n=l 

= / (/o2A/(i+i)i - f(i+i)iAfo2)dn 

JSG 



/ , 



f02flldfi- 
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This shows that our definition is consistent with [SU]. It is easy to see that a_i = 
2,6_i = l. 

We note here some typos from [SU] : 

(i) in (5.4) the coefficient || should be ||; 

(ii) in the first fine of (5.7) the coefficients 2 of aj-i-£ and hj-i-i should be 
deleted. 

Now let Pj, qj be defined by 

Pj = ^^fjk{Fiqk) i 7^ k 

^3 = fjk {Fiqd for ^5 3 1 ^ distinct. 

(Note that we are using the same symbol qj for two different things, but it should 
be clear from context which is which.) 
Then (5.7) of [SU] rearranged says 

3 

^(aj_/_i + bj_i_i){2pi + qi) + = 

1=0 
3 

J^(2oj_/_i - hj-i-i){pi - qi) + hj-i = 
/=o 

If we set 



A = 



/a_i 























do 


a_ 


1 











bo 


b-i 







ai 


Go 




a_i 







B = 


bi 


bo 


b-i 







ai 




ao 


a_i 






b2 


bi 


bo 


b-i ■•• 


V ; 














\ \ 








■J 




/Po 









\ 




(qo 







\ 






Pi 


Po 











qi 


qo 








p = 


P2 


Pi 


Po 







Q = 


q2 


qi qo 











P3 


P2 


Pi 


Po 






qs 


q2 qi 


qo 








\ : 












V : 






••J 





Then in matrix notation this becomes 



(2.13) {A + B){2P + Q) + B ^0, 



{2A- B){P -Q)+B ^ 0. 
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Now for j >0, 



(2.14) 



1=0 

j 



1=0 



SO taking normal derivatives at go) we have 



-1 + 2^aj_i6i_i = dnPjiiqo) = 



2j2/3j-lbl-i=dnPj2iqo) = 



1=0 



In matrix notation this is 



1, ifj=0; 
0, otherwise. 



I.e. 



(2.15) 



2aB + A = 0, 2(3B = /, 



Substituting (2.15) into (2.13), we get 

2P + Q = -{A + B)-^B = -[-l/3-i(2a - /)]-M^/3-'] = (2« - /) 
P-Q = -{2A - B)-'B = -[-i/3-i(4a + 1)]-^^^ = (4« + 



so 



{2a - /)(2P + Q) = / = (4a + I){P - Q). 



Expanding we get 



AaP + 2aQ - 2P - Q = AaP - AaQ + P - Q, 



I.e. 



(2.16) P = 2aQ, and Q = (Aa + I)-\2a - 

Now evaluate (2.14) at Fq^i, noting that 

Pji{Foqi) = 5-'Pji{n) = 
Pj2{Foqi) = h-^PjiiQi) = l^-'f^j 
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by (2.4), (2.5) and 

^2 floiFoqi) = fii{Foqi) = 5"^^ 



by the definitions of pi^s and g;'s. The result is 



1=0 

J 



1=0 



so 



= ^'j + XI ^(^j-i{Pi + (li) and 

1=0 



In matrix notation these read 



From (2.14) we see that 



a = P + t(q;)(P + (3) and 

^/3 = T(/3)(p + g). 

o 



a = [2q; + T{a){2a + I)\Q and 

^/3 = t(/?)(2« + /)Q 
5 



hence 



t(q;) = 4q!'^ — 3q; and 

|/3(2a - /)(4q! + /) = T(/3)(2a + /), 
5 

from which (2.9) and (2.11) foUow. 
Finally 

j 

= Ylj-i{fii- fn) 

1=0 
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and so by (2.17) we have 

1=0 

i.e. 



1 ^ 



1=0 

or in matrix notation 

^7 = t(7)(P-Q). 

5 

Thus t(7) = I (4a + 1)7 from which (2.10) foUows. 

The values of ao, Po and 70 are easy to check. Then (2.12) follows from (2.9) 
and (2.10) since aj and ckj-i satifsy the same recursion relation. Q.E.D. 



Theorem 2.4. For all j >0 we have 



(2.18) Pl^\x) + Pff{x) + Pff{x) = 



and 

(2.19) Pi!\x) = 3(P^,^,(x) - P(^i)i(x)). 

Proof : We prove (2.18) by induction. For j = the left side is a harmonic 
function that vanishes on the boundary (because of the skew-symmetric of each 
term). Such a function must be zero. For the induction step, assume it is true for 
j -1. Then 

Afp(0) + p(l) + p(2)) _ p(0) , p(l) , p(2)_Q 
^K^j3 + + ) - -^0-1)3 + -^0-1)3 + ^j3 - U 

by the induction hypothesis. Once again the left side is a harmonic function, and 
it vanishes on the boundary by skew symmetry. 
To prove (2.19) we use 

(2.20) P^l^ = j2^j-i{fn-fi2). 

e=o 

On the other hand, we have 

i+i 



P(j+i)i = /(i+i)2 + ^(Xj-£+i{fm + In) 



>(2) 

£=0 

j+i 



= /o+i)i + ^oij-e+i{f£o + fa) 
e=o 



16 J. NEEDLEMAN, R. S. STRICHARTZ, A. TEPLYAEV AND P-L. YUNG 



SO that 

i+i 



~ - /(i+i)2 - /(i+i)i + - /e) 

3 

= ^^oij-i+i{fn - fa) 



£=0 



since ao = 1. The result follows from (2.12). Q.E.D. 

The dihedral-3 symmetry group -D3 of SG consists of reflections po, pi, p2, where 
Pj preserves qj and permutes the other two boundary points, and rotations /, 
R2 = (Ri)^ where RiQj = Qj+i (cyclic notation). 

Theorem 2.5. Any polynomial P satisfies the identity 

(2.21) P{x) + P{Rix) + P{R2x) = P{pox) + P{pix) + P{p2x), 
and more generally the local versions 

(2.22) P{xo) + P{x^) + P{x2) = P{yi) + P{y2) + P{yz) 
for any sextuplet of points such that 

^2 I ^0 = ^1 = KRix, X2 = F^R2X, 

I yo = RwPox, yi = F^pix, y2F^p2X 

for some x E SG and some word w. 

Proof : The local version follows from (2.21) because P oFyj is also a polynomial. 
To prove (2.21) it suffices to show it holds for all monomials. Now we claim (2.21) is 
trivially true for any function that is symmetric with respect to one of the reflections 
Pj. Say P{x) = P{pox) for aU x. Then P{Rix) = P{pix) and P{R2x) = P{p2x) 

because po-Ri = Pi and P0-R2 = P2- In particular, (2.21) holds for all P^^^ and Pj2- 

It follows from (2.19) that it also holds for Pf^. Q.E.D. 

The same result holds for uniform limits of polynomials; in particular, the con- 
vergent power series discussed in the next section. Note that Kigami [Ki2] Theorem 
4.3.6 has characterized the space of limits of polynomials by the condition of 
orthogonality to all joint Dirichlet and Neumann eigenfunctions. It is not hard to 
see that (2.22) implies the orthogonality to some of these eigenfunctions (those of 
the A*^^^-type in [DSV]), but not others. On the other hand, it is not clear how 
these orthogonality conditions imply (2.22). 
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Corollary 2.6. Any polynomial P satisfies 

(2.24) drPigo) + drPiqi) + dTP{q2) = 0, 

and more generally the sum of tangential derivtives at the boundary points of any 

cell must vanish. 

Proof : Taking x = Fq^qi in (2.21), we find 

(2.25) {P{FS'q^)-P{F^q2)) + {P{Fr'l2)-P{Frqo)) + {P{F^qo)-P{F^ql)) = 

because RiF^qi = ^^^2, R2F^qi = F^qo, poF^qi = F^q2, piF^qi = F^qi, 
P2F^qi = F^qQ. Multiplying (2.25) by 5"^ and taking the limit as m — > oo yields 
(2.24). The local form follows as before. Q.E.D. 

Remark : As we observed in the proof of Theorem 2.5, any polynomial may be 
written as a sum of three polynomials, each symmetric with respect to one of the 
refiections pj, P = P^^^ + P^^^ + P^'^\ It is easy to see that one way to do this 
explicitly is to take 

(2.26) P^^\x) = ^-{P{x) + P{p^x)) - \{P{pox) + P{pix) + P{p2x)). 

We consider next estimates for the size of aj , fij , 7^ . We show that aj has rapid 
decay, which we believe is fairly sharp. This gives the same decay rate for 7^-. 

Theorem 2.7. There exists a constant c such that 

(2.27) 0<aj <c(j!)-'°s^/i°s2 j 

Proof : It is clear from (2.9) and the initial conditions that the a.j are positive. 
Let aj = (j!)i°s^/i°82a^.. We need to show that the aj are bounded, which we do 
by induction. If 5^ < c for £ < j, then (2.9) implies 



i-l / -\ log 5/ log 2 

a,.cv-.i:(3 



It is well known that 

j 

E 



so by Stirling's formula and routine arguments we have 

i-l / A log 5/ log 2 



<M5^(i)- 



for all j > 2 for a small constant M, so aj < c^5M{j) ^1'^ . It is easy to choose c 
and jo so that 5^ < c for £ < jo and c < {joy/'^/5M. Q.E.D. 

Table 2.1 presents numerical computations of aj and Pj. 
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J 


aj 




I3j 






. ,^ log(5) 





1. 




-.5000000000 




- 5000000000 

■ \j W W \j \j \J \j \J \J w 


1. 


1 


.1666666667 




-.04444444444 


6 025427867 


1 333333333 


2 


.005555555556 


-.001008230453 


-18 53107571 

-L \^ • \_/ \J -1- V/ 1 \_/ 1 -J- 


1 111111111 

-t-*lllllllll 


3 


.00006172839506 


-.8554950809 10 


_5 


21 31713060 


2 025658338 


4 


.3318730917 10 


-6 


-.3853047646 10 


_7 


-13.01625411 


2.178127244 


5 


.1021147975 10 


-8 


—.9848282711 10 


— 10 


4 510374011 


2 250339083 


6 


.2007235906 10 


-11 


—.1933836698 10 


-12 


-1 200721414 


2 268082964 


7 


.2713115918 10- 


-14 


—.7720311754 10 


-16 


06498718216 


2 248411184 


8 


.2656437390 10" 


-17 


-.1187366658 10" 


— 17 


- 1355027558 


2 201440598 


9 


.195916520110" 


20 


.7232200062 10" 


20 


- 1118933095 


2 134277683 


10 


.1122370097 10" 


-23 


-.5436238235 10" 


-22 


- 1140256558 


2.052740417 


11 


.5120236416 10" 


-27 


.4004514705 10" 


24 


- 1138739539 


1.961629028 


1 9 


.1898528071 10" 


-30 


-.2954013973 10" 


-26 




J- .OUt: I ZiUtit:!- 


13 


.5820142006 10" 


-34 


.2178916451 10" 


28 


-.1138822148 


1.764891613 


14 


.1496625756 10" 


-37 


-.1607201123 10" 


-30 


-.1138822304 


1.664234594 


15 


.3268360869 10" 


-41 


.1185495242 10" 


32 


-.1138822298 


1.564302197 


16 


.6126918156 10" 


-45 


-.8744387717 10" 


-35 


-.1138822298 


1.466232140 


17 


.9952451630 10" 


-49 


.6449989323 10" 


37 


-.1138822298 


1.370864839 


18 


.1412543698 10" 


-52 


-.4757607235 10" 


-39 


-.1138822298 


1.278818576 


19 


.1764707126 10" 


-56 


.3509281252 10" 


41 


-.1138822298 


1.190538877 


20 


.1953558627 10" 


-60 


-.2588497599 10" 


-43 


-.1138822298 


1.106332006 



Table 2.1. 

It appears that 8-'(j!)'°^^/^°^^aj remains bounded (8 is by no means the best 
constant, and perhaps it could be replaced by an arbitrary positive number). 
It also appears that {—\'2)^iij converges to the constant —.1138822298, where 
A2 = 135.572126995788 ... is the second nonzero Neumann eigenvalue. It is easy 
to see that A2 is the largest value for which such an estimate could hold, because 

00 

^/?j(-A2)^ diverges. 

j=0 

Indeed, if we did not have divergence then 

00 

would be a solution to the eigenvalue equation —Au = X2U satisfying dnU^qo) = 1. 
But, since A2 is not a Dirichlet eigenvalue, the space of eigenfunctions has dimension 
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three, whereas the multiphcity of the A2-Neumann eigenspace is also three, so every 
eigenfunction automatically satisfies dnu{qo) = 0. 

We note that the computation of Pj, carried out using the recursion relation 
(2.11), was done using exact rational arithmetic (the reported values are reported 
as decimal approximations, of course). This is significant because this solution of 
(2.11) is highly unstable. For example, if we take Po = ^ and (3i = .044444444 or 
.04444445 (the correct value being 2/45) and then use (2.11) for j > 2, we find the 
ratio f3j/(3j+i approaching -84.0799 . . . (this is -5Af , where Af = 16.815999 . . . 
is the first Dirichlet eigenvalue). In Section 6 we will give an explanation for this 
phenomenon. 

Next we will establish estimates for ||Pjfc||oo- To do this we will study the oper- 
ator 

(2.28) Af{x) = Gf{x) - {dniGf){qo))Po2 

where Gf{x) = j G{x,y)f{y)dii{y) is the Green's operator, satisfying —A.Gf = f 
and Gf{qi) = 0, z = 0,1,2. Note that A is a compact linear operator, but is 
not self-adjoint. Thus the spectrum of A consists of isolated eigenvalues of finite 
multiplicity, and zero. Note that we have 

(2.29) -AAf = /, Af (go) = and d^Af (go) = 0. 
In particular, this implies 

(2.30) AP.fc = for A; =1,2. 

Write Aq for the restriction of A to the i?o-symmetic functions, where Rq is the 
refiection preserving go. 

Lemma 2.8. (a) f is an eigenfunction of Aq {Aof — A/) if and only if f is a 
symmetric -eigenfunction of A satisfying /(go) = dnfiqo) = 0. (b) f is an 
eigenfunction of Aq if and only if f is a symmetric -Neumann eigenfunction of 
A satisfying f{qo) — 0- (c) The Jordan block of Aq associated to any eigenvalue is 
diagonal. 

Proof : (a) By (2.29), any eigenfunction of ^ is a A~ ^-eigenfunction of A satis- 
fying f{go) = dnf {go) = 0. For the converse, let v = Af — A/. Then 

= AAf - XAf = A{Gf - dn{Gf)P2) +f = -f + f^O 

so V is harmonic. But v is symmetric with v{go) = dnv{go) = 0, and this implies 
v = 0. 

(b) The only new assertion here is that / in part (a) also satisfies dnf{Qi) = 
dnf{q2) = 0. This requires a rather detailed knowledge of the description of 
eigenfunctions of A by spectral decimination. First we observe that if |A~^| is 
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small enough (less than the first Dirichlet eigenvalue), then a symmetric A~^- 
eigenfunction is uniquely determined by f{qo) and dnf (go)- This implies that / 
vanishes identically on a cell Fq{SG) for n large enough. But an eigenfunction 
can vanish on a cell only if the space of eigenf unctions has dimension greater than 
three, and that happens only if is a joint Dirichlet-Neumann eigenvalue. That 
means its restriction to the graph for some value of m is either a 5-eigenfunction 
or a 6-eigenfunction. In the 6-eigenfunction case there is nothing to prove, since 
all eigenfunctions are Neumann eigenfunctions. In the 5-eigenfunction case this is 
not true, but the Neumann eigenftmctions have codimension two in the space of 
all eigenfunctions. When we impose the -Ro^symmetry condition the codimension 
drops to one. We know exactly what this one function looks like (see Figure 2.2 
for the case m = 2). In particular, it does not vanish identically in any small cell 
Fq^{SG). Since / does (and so do all symmetric joint Dirichlet-Neumann eigen- 
functions), it follows that / must be Neumann eigenfunction (in the 5-eigenfunction 
case it is also a Dirichlet eigenfunction, but not necessarily in the 6-eigenfunction 
case). 



Figure 2.2. 

(c) Suppose A is an eigenvalue of Aq, and {Aq — X)^g = 0. Then A~^ is a 
Neumann eigenvalue of A, and (A + X~^)^g = 0. Also g is symmetric and satisfies 
9iQo) = dng{qo) = 0. By similar reasoning as before, ^ is a Neumann eigenfunction 
of A, hence the Jordan block associated with A is diagonal. Q.E.D. 

Theorem 2.9. (a) For any r < oo there exists Cr such that 

(2.31) ||^':?i||oo < Crr-\ 
or more precisely 

(2.32) lim - log ||P,i||oo = -oo. 

i^oo j 

(b) There exists c such that 

(2.33) ||Pi2||oo <cA-^ 
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and 

(2.34) lim (-A2)^Pi2 = ^ 

j—^oo 

where (p is a \2-Neumann eigenf unction of A which is Ro-symmetric and vanishes 
on Fq{SG) (a multiple of the eigenfunction shown in Figure 2.3 on Ti), the limit 
existing uniformly and in energy. 

Proof : (a) Consider the norm 

(2.35) 11/11 = (11/11^ + ^(/, 

and define Ci and £2 as the closures in this norm of the spans of {Pji} and {Pj2}, 
respectively. By (2.30), Aq preserves both spaces. Denote by Ai and A2 the 
restriction of ^0 to jCi and £2- We claim cr(^i) = {0}. Indeed, otherwise Ai would 
have to have a nonzero eigenvalue A because Ai is compact. Since this would also 
be an eigenvalue of Aq, by Lemma 2.8 would have to be a Neumann eigenvalue 
of A. So A > 0, and we may choose it to be the largest eigenvalue of Ai. Then 
X~^A\ converges to a projection (not necessarily orthogonal) Bx onto the finite 
dimensional A-eigenspace of Ai. Note that BxPqi cannot be the zero function, 
because that would imply BxPji = for all j, contradicting the fact that Bx is 
nonzero. But then A-[Pqi = Pji would converge to a nonzero eigenfunction 
of Ai. By Theorem 2.7 this eigenfunction would vanish at qi and q2, and of course 
it vanishes at qo, since Pji does for j > 1. So it would have to be a joint Dirichlet- 
Neumann eigenfunction of A. But Theorem 4.3.6 of [Ki2] asserts that all Pjk are 
orthogonal to all joint Dirichlet-Neumann eigenfunctions. 

Thus we have shown that o-{Ai) = {0}, so the spectral radius of Ai is zero, 

lim \\A{\\^/^ = 0. 

Applying this to Pqi we obtain (2.32) (the norm (2.35) dominates the L°° norm), 
which implies (2.31). 

(b) The result of Kigami used above moreover says that £ = £1 © £2 con- 
tains all i?o-symmetric Neumann eigenfunctions of A that are othogonal to all 
joint Dirichlet-Neumann eigenfunctions (note that Kigami uses the norm rather 
than (2.35), but the same argument applies). In particular, it contains the A2- 
eigenfunction shown in Figure 2.3 (this is a Neumann eigenfunction, so it is or- 
thogonal to all Neumann eigenfunctions with different eigenvalues, and there are 
no joint Dirichlet-Neumann eigenfunctions with the same eigenvalue). By Lemma 
2.8 and the explicit description of Neumann eigenfunctions, A^"^ is the largest eigen- 
value of Aq, and ip spans this multiplicity one eigenspace. Thus, as before, X2A^ 
converges to a one-dimensional projection operator B^-i, and S^-iPoi = 0. That 
means B^-iPo2 7^ 0, for otherwise B^-i — 0. So 

lim (-A2)^P,-2 = lim X^A^ P02 = B.-1P02 

j—fOO j^OO 2 
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which is (2.34). This imphes (2.33). 



Q.E.D. 



The estimate (2.33) is sharp, but (2.32) falls short of what we would have if we 
knew ll-Pjilloo = ctj, in view of (2.27). One approach to establish this would be to 
prove the following conjecture: 



We have numerical evidence for this conjecture for moderate values of j. To 
show that (2.36) implies ||-Pji||oo = ctj is easy using the following well-known fact 
(we provide a proof since it does not appear explicitly in the literature) . 

Proposition 2.11. If u E dom A, Att(xo) > and xq is not a boundary point, 

then u does not achieve its maximum value at xq. 

Proof : If Xq is a vertex in 14 the result follows immediately from the pointwise 
definition of /S.u{xq). If not, then we can find a cell F^K such that xq is in the 
interior of F^K and Ait > on F^jK. Let v = uo F^. Then Av > 0, and we have 



where G is the Dirichlct Green's function and h{x) is the harmonic function with 
the same boundary values as v{x). Since the Green's function is positive in the 
interior, we have v{x) < h{x) in the interior. Since h attains its maximum on the 
boundary, it follows that v cannot attain its maximum in the interior, so u{xo) is 



Figure 2.3 



Conjecture 2.10. For all x ^ qq and all j, 



(2.36) 



P,i(a;)>0. 




not a maximum. 



Q.E.D. 



Next we study the normal and tangential derivatives of monomials at boundary 
points. 
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Theorem 2.12. We have initial values uq = 0, to — —1/2, and recursion relations 

j-i 

(2.37) nj = ^^^aj + 2j2ne(3j-£ for j>l, 

e=o 



i-i 

(2.38) tj^(3j-6j2aj+i-ete for j>l. 



Moreover, we have 



(2.39) dnPj2{qi) = dnPj2{q2) = | 



I - ao if j = 
-oij if j > 1 



(2.40) dnPjsiqi) = -dnPoM = 3n,+i 



(2.41) dTPji{qi) ^ -dTPji{q2) 



(2.42) dTPjz{qi) = -dTP^^iq^) = | 



-e if 3 = 1 
^/ jVI 



-\ ^f J = 
^/ i>l 



Proof : As in the proof of Theorem 2.3 we introduce matrices n, n and t, where 
rij = dnPj2{qi)- When we evaluate the normal derivatives on both sides of (2.14) 
at qi, we see that 



rij = bj-i + aj-i{ai^i + 6/_i) for aU j, 

1=0 

or in matrix notations 

n = B + a{A + B). 

Using (2.15) this yields 

(2.43) n = -I3~\I + 2a){I - a) ^ 7/5"^(2/ - T{a) - a) 

which implies (2.37). 
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By the same reasoning 

j 

Uj = ^f3j-i{ai-i + bi-i) for all j. 

1=0 

Then 

n = P{A + B) 

and hence by (2.15) we obtain 

(2.44) n=\l-a, 

which implies (2.39). 

Finally, the same reasoning shows 



tj = J2Pj-iTi for all j, 



where T/ = drfaiqi)- Now Pjs = Yl\=o lj-i{fii - fn), so taking tangential deriva- 
tives at qo we get 



In matrix notations these become 



1, ifj = 0; 
0, otherwise. 



t = (3T 

Together we have 

(2.45) /5 = 27t = 6(j(a)t, 

where the last equality follows form (2.12). 

This proves (2.38). The initial values of no, and to arc easy to check. 

Note that the skew-symmetry implies dTPjsiqi) = dTPj3{q2), so (2.2) implies 
drPjsiqo) + ^drPjsiqi) = 0, which yields (2.42). Then (2.41) follows from (2.19) 
and (2.42), and similarly (2.19) implies (2.40). Q.E.D. 



CALCULUS ON THE SIERPINSKI GASKET 



25 



Theorem 2.13. For any r < oo there exists such that, for all j >l, 

(2.46) \nj\ < Crr-^. 
Also 

(2.47) \tj\<cX2^. 



Proof : From the Gauss-Green formula we have 

2 

/ Audu = ^dnu{qi). 

We apply this tou = PjJ^ , noting that a^Pj?^ (go) = and d^P^^^^ (gi) = a^Pj?^ (ga) = 
Uj. It follows that 

(2.48) ""'^Ij ^(i-i)!^^' 

and (2.46) follows from (2.31). 

Similarly, (2.47) will follow from (2.33) and the estimate 

(2.49) \dTu{qi)\ < cdl^lloo 

In [S3] it is shown that dTu{qi) exists if u e domA and Au satisfies a Holder 
condition, and (2.49) is just a quantitative version of this fact. For the convenience 
of the reader we outline the argument. For simplicity take i = 0. Let Qm (see 
Figure 2.4 for m — 2) denote the level m piecewise harmonic function satisfying 
gmiqo) = and gm{F^qi) = 3^= and gm{F^q2) = -3^= for aU k < m. Then 



Figure 2.4 
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(2.50) J gmAudf, = ^5"'{u{F^q^) - u{F^q2)) - h{u{q^) - u{q2)) 



by the Gauss-Green formula, since the sum of the normal derivatives of at Fq^qi 
is (14/3)5"^ (there are no terms involving normal derivatives of u at Fo^qi because 
u satisfies matching conditions). Let ui = Au. Note that gm is odd, so only the 
odd part of ui contributes to the integral in (2.50). So (2.49) will follow from (2.50) 
and the estimate 



(2.52) J gm{ui - ui o RQ)diJ, < c(||ui||oo + ||Aui||oo)- 



But (2.52) is routine, because on the cells F^Fi{SG) and F^F2{SG) {0 < k < m) 
of measure S~''~^, the function g^ is of size S'^, and ui —uio Rq can be estimated 
by (|)'=||A«i||oo. Q.E.D. 



In Table 5.2 we display the results of solving the recursion relations for nj and 
tj. The data suggests that (— A2)"'tj converges, in fact quite a bit faster than for 
(3j, and lim I3j/tj+i — 9. Moreover n,- is always positive and satisfies 



(2.53) rij < cjaj. 



If Conjecture 2.10 holds, then ||-P(j-i)i ||oo = cxj-i so (2.48) implies nj < ^aj-i, 
which is only slightly weaker than (2.53). 
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.23573795 


10- 


-47 


-.97159864 10-3^ 


1.3933188 


1.7154736 


9.0000000 


18 


.34893132 


10- 


-51 


.71666548 10-3^ 


1.3723521 


1.7154736 


9.0000000 


19 


.45359082 


10" 


-55 


-.52862303 10'^^ 


1.3528138 


1.7154736 


9.0000000 


20 


.52141937 


10" 


-59 


.38992014 10-42 


1.3345373 


1.7154737 


9.0000000 



Table 2.2. 



We also have found that the recursion relation for nj is unstable, and any slight 
perturbation produces a decay rate 0((Af )-^), which is even slower than the decay 
rate for Pj and tj. Also a slight perturbation of the tj recursion relation produces 
a decay rate of ©((Af)-^). We will explain this in Section 6. 

Next we describe the change of basis formula to pass between {Pjj^} for different 
values of ^, an immediate consequence of Theorem 2.12. 



Corollary 2.14. We have 



(2.54) 



/ -Ji 

\^j3 



\^k3 / 
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for matrices Mj given by 



(2.55) 



Similarly 

(2.56) 

for 



(2.57) 







nj 


o\ 










\3aj+i 


-aj 
3nj+i 


^) 


for 


J >2 






/ ai 


ni i] 






/ ao no 


\ 


Ml = 


Pi 


—ai tl 




Mo = 


(3q \ — ao 


to 




\Sa2 


3n2 J 






\ ScKi 3ni 


-u 



j-k 



k=0 



p(^-l) 
\-^fe3 / 





/ 




nj 

























-aj 








for 


i >2 








I 




-3nj+i 





) 
















/ 




ni -\\ 










f 




no 


\ 


Ml = 






— tti — tl 




Mo 






/3o 








V 


-3q;2 


-3n2 ) 












— 3q;i 


— 3ni 


-u 



§3. Power series. 

A formal power series about qi is an expression of the form 

3 oo 

(3.1) E J2^jkP^k(^)- 

k=l j=0 

We call {cjk} the coefficients, and we seek growth conditions on the coefficients 
that will make (3.1) converge nicely. 

Theorem 3.1. If the coefficients satisfy 

(3.2) |cji| and \cj2\ = 0{{jiy) for some r < log 5/ log 2, 
and 

(3.3) \cj2\ = 0{R^) for some R < X2 
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then (3.1) converges uniformly and absolutely to a function u e dom{A°°), and 
(3.1) may he "differentiated term-by-term", 

3 oo 

(3.4) ^M^) = Y: E^.^^O-n).(^)- 

fc=l j=n 

Moreover, the coefficients are given by the infinite jet of u at qg: 

{Cji = A^u{qe) 
Cj2 = dnA^u{qe) 
Cj3 = dTA^u{qi). 

Proof : The estimates in Theorem 2.9 conspire with the growth rates (3.2) and 
(3.3) to make (3.1) converge uniformly and absolutely. Call the limit u. Note that 
the right side (3.4) is also a formal power series, in fact 

3 oo 
k=l j=0 

whose coefficients also satisfy the growth rate conditions (3.2) and (3.4). So the 
right side of (3.4) converges uniformly and absolutely. By terminating the sums at 
j = N and letting AT — > oo we obtain the equality in (3.4) by a routine argument 
using the Green's function [Ki2]. 

It suffices to prove the jet formulas (3.5) when j = in view of (3.4), and for this 
it suffices to show that if cqi = co2 = cqs = then u{q£) = dnu{q() = dru^qi) = 0. 
Of course u{q() = directly from (3.1). For simplicity put £ = 0. Then (since 
<qo) = 0) 

a,«(go) = - lim (-) {u{F^q^) + u{F^q2)). 

But we have 

oo „ 

(3.6) u{F^x) = J2cji5-^'Pji{x) + c,2{-5-^y Pj2{x) + c,35-^(^+'^ P,s{x). 

Using the estimates for the coefficients and monomials we see that 

(3.7) u{F^x) = 0(5—), 

and this suffices to prove dnu{qo) = 0. This by itself does not suffice for the tan- 
gential derivative, which has a factor of 5"^. However, for the tangential derivative 
we can restrict attention to the skew-symmetric part 

^ oo 

(3.8) u{x) = -{u{x) - u{pqx)) ^^Cj^Pj^ix), 
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SO the analog of (3.6) shows 

(3.9) u{F^x) = 0(5-2^), 

which implies dTu{qo) = 0. Q.E.D. 

As a corollary of the proof we can characterize rates of vanishing of power series. 

Definition 3.2 : A function / is said to vanish to order r (any positive real) at q£ 
provided 

(3.10) ll/oFf |U = 0(5— ^). 

If (3.10) holds for all r then we say / vanishes to infinite order at qg. 

Corollary 3.3. If u is represented by a power series (3.1) with coefficients sat- 
isfying growth conditions (3.2) and (3.3), then u vanishes to order N (a positive 
integer) at qi if and only if cjk — for all j < N. In that case A^u vanishes to 
order N — i for all £ < N. Moreover, the odd part u vanishes to order N + 1. In 
particular, if u is not identically zero then it cannot vanish to infinite order. 

Next we consider rearrangement of power series, moving from one boundary 
point to another. It turns out that we need to make stronger assumptions on 
the coefficients, requiring Cji and c^s to satisfy the same exponential growth rate 
as Cj2- 

Theorem 3.4. Suppose the coefficients of a power series (3.1) about one boundary 
point qi satisfy 

(3.11) \cjk\ = 0(W) for some R < X2, k = 1, 2, 3. 

Then the function may also be represented by power series about the other boundary 
points with coefficients also satisfying (3.11). More precisely, the coefficients at g^+i 
are given by 

00 

(3.12) (C^.,1 C^v2 C^vg) - ^(C(j+jv)i C(j+jv)2 C^j+j,)s)Mj 

and similarly at q£-i with Mj replaced by Mj (see (2.55) and (2.57)). 

Proof : The key observation is that the right side of (3.12) converges absolutely 
and the new coefficients again satisfy (3.11) (in fact with the same value of R) 
because the entries in Mj are ©(A^"') by Theorem 2.13. Of course (3.11) is ex- 
actly what we get if we substitute (2.54) into (3.1) and interchange the order of 
summation, which is easily justified using the estimates of Theorem 2.9. Q.E.D. 
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Note that we could not allow slower growth rates like (3.2) for the Cji and 
Cjs coefficients and still rearrange, because the second column of Mj has positive 
entries. In §6 we will present an example to show that rearrangement fails when 
Cji = 0(A2). However, condition (3.11) is not sharp. We could replace it by 

oo 

(3.13) 5Z-^2'^|CjA;| < OO, 

j=0 

and the rearranged coefficients would satisfy the same growth condition. However, 
not all subsequent results would be valid under this hypothesis. 

Definition 3.5 : An entire analytic function is a function given by a power series 
(3.1) with coefficients satisfying (3.11). 

We can also consider local power series expansions on any cell Fyj{SG) with 
respect to a boundary point F^jQi of the cell, namely 

oo „ 

(3.14) Yl {^-""'^^nPi^iK'^) + {f-'yc,,P^2\F-'x) 

+ 5-0-+^)-c,3P,^?(F-M) 

where m = \w\. 

Theorem 3.6. An entire analytic function has a local power series expansion 

(3.14) for any w and £ with coefficients satisfying (3.11). Conversely, suppose u{x) 
is a function defined on Fyj{SG) given by a local power series expansion (3.14) with 
coefficients satisfying (3.11). Then u has a unique extension to an entire analytic 
function. 

Proof : Suppose first that m = 1, say w — iff). If £ = then the local and global 
power series are identical, with identical coefficients. Moreover, u o is an entire 
analytic function with coefficients satisfying (3.11) (in fact with R < A2/5). The 
rearrangement for u o about qi and q2 guaranteed by Theorem 3.4 gives the 
local power series of u in ^0(5*^) about -FqQ'i and -FoQ'25 with the same coefficient 
estimates. We may then iterate this argument to get local power series about any 
boundary point in any cell. 

Conversely, suppose u is given in F^{SG) by a local power series about Fy^qg, 
with coefficients satisfying (3.11). Write w = iw' , Wm) with \w'\ = m — 1. If Wm 7^ £ 
then use Theorem 3.4 to rearrange the power series of w o F^ about qw^- So we 
end up with a local power series of u about F^'Fgqg in the cell FyjiFi[SG). But 
Fyj'Fiqi = Fw'qe and the power series makes sense in the cell Fyj'{SG). Use this 
power scries to extend the definition of u. By iterating the argument, we obtain the 
desired extension. Note that the estimates (3.11) on the coefficients are reproduced 
in each extension or rearrangement step. It is clear that the extension is unique 
because the rearranged coefficients are determined by (3.12). Q.E.D. 
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By the same reasoning, if a local power series has coefficients satisfying 

(3.15) Cjk = 0(W) for some R < 5^°A2, 

then the function can be also represented by a power series on a level uiq cell. 
One might hope that this "analytic continuation" might extend somewhat beyond 
the cell, with the domain of analyticity growing as R decreases toward 5'^°~^X2. 
However, the experimental evidence we have seen does not support this at all. On 
the contrary, we will see in §6 that there are power series (3.1) with coefficients 
0(A2) where we have divergence outside Fe{SG). We might describe this as a 
"quantized radius of convergence." Of course, this does not rule out a different 
type of behavior for special classes of power series. 

Theorem 3.7. An entire analytic function satisfies the estimate 

(3.16) IIA'^mIIoo = OiR"") for some R < Ag. 



Proof: We have 



3 oo 

{j-n)k 

k=l j=n 



so 



3 oo 



iiA"«iu < E^'ii^u-n)feiioo < Mj2R'>^r' = 0{Rn 

k=l j=n j=n 

for R in (3.11). Q.E.D. 

The condition (3.16) obviously implies the same estimate in norm: 

(3.17) ||A"w||2 = OiR"") for some R < A2. 

But conversely, (3.17) implies (3.16), because ||/||oo < c(||/||2 + IIA/II2). The 
estimate (3.17) is technically more convenient, since we can compute norms 
exactly from eigenfunction expansions. 

It follows immediately from the definition that an eigenfunction of A is an entire 
analytic function if and only if the eigenvalue satisfies |A| < A2. Theorem 3.7 shows 
us that many other functions that we might believe to be entire analytic functions 
are not. Indeed, suppose u is represented by a Dirichlet (or Neumann) eigenfunction 
expansion 

00 

(3.18) u{x) = ^ak(pk{x) 

k=l 
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where {^Pk} is an orthonormal basis of Dirichlet (or Neumann) eigenfunctions. If 
the coefficients are rapidly decreasing, 

(3.19) ak = Oik-"") for all n, 
then we may differentiate term-by-term, 

oo 

(3.20) AX^)=^(Af)'^afe(^fe(x). 

k=l 

It follows that 

°° 1/2 

(3.21) ||A-«||.= (5^(Af)^-|a,r) . 

fc=i 

If (3.18) is non-trivial in the sense that an infinite number of coefficients are non- 
zero, then not only does (3.17) fail to hold, but the estimate cannot hold for any 
finite R. So u cannot be represented by a local power series with (3.14) holding on 
any cell. In particular this applies to the heat kernel. 

This observation stands in striking contrast to the situation on the unit interval, 
where analyticity properties of a function may be characterized by decay properties 
of the coefficients of its Fourier series expansion. 

§4. Characterization of analytic functions. 

The main purpose of this section is to prove the following theorem. 

Theorem 4.1. u is an entire analytic function if and only if u & dom{A°°) and 
(3.16) (or equivalently (3.17)) holds. 

We first consider the case when u is even with respect to po- In that case we 
would like a Taylor expansion with remainder about qo, 

(4.1) u{x) =Tku{x)+Rk{x) 
for 

k-i 

(4.2) Tku{x) = J2 ^'u{qo)Pji{x) + {dnA^u{qo))Pj2{x) 

3=0 

and Rk{x) the remainder term. While we can use (4.1) to define the remainder, to 
be useful we need some explicit expression for it. We are only able to do this for 
X ^qi (or q2). 
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Lemma 4.2. Let be a function in Hk that is even with respect to po satisfying 
(4.3) A^Vk{qi)=Oforj<k-l 



(4.4) dnA^Vkiqi) 
Then 



for j <k-2 
-\ forj = k-l. 



(4.5) Rk{qi) = RM = / VkA'^udii 

JSG 

for even functions u e dom{A^). 

Proof : Note that A^u = A^{u - T^u) = A^Rk. We apply the Gauss-Green 
formula k times to obtain 



j VkA^udii. = j 



VkA'^ud/j, = I VkA^Rkd^ 

k-i 

= 2Y,{^'vk{qi)dnA''-^-'Rk{qi)-dnA^Vk{qi)A>'-^-'Rk{qi)) 
j=o 

since A^-^-^Rk{qo) = dnA^-^-^{qQ) = 0. By (4.3) and (4.4) all terms vanish 
except when j = k — 1 and we obtain exactly Rk{qi)- Q.E.D. 

Lemma 4.3. The function 

k-i 

(4.6) Vk = Y^i-Pk-l-iP^i^ + ak-i-iPg^) 

satisfies the conditions of Lemma 4-2. 

Proof : Clearly Vk G 7ik and is even. Since 

k-j-i 

A^Vk = -Pk-j-i-iPn + ak-j-i-iPa 

£=0 

we obtain 

k-j-i 

^^Vk{qi) = ^ {-Pk-j-i-£Oie + ak-j-i-iPe) = 
which is (4.3). Similarly 

k-j-l k-j-1 ^ 

d„AHk{qi)= ^ Pk-j-i-im- ak-j-i-iai^ + -ak-j-1 
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by (2.35). When j = k-l this is just 



For j < k — 2 we have 

Pk-j-i-eni = 



2^ Pk-j-i-en£ = -[^ )ak-j-i 



by (2.33), and 



1=0 

by (2.9) (this uses k-j-l>l). Thus dnA^Vkiqi) = 0, proving (4.4). Q.E.D. 

Lemma 4.4. If u is an even function in dom{A^) satisfying (3.16), and u is the 
entire analytic function whose expansion about qo has coefficients cji — A^u{qo), 
Cj2 = dnA^u{qo) and Cjs ~ 0, then u{qi) — u{qi) and u{q2) — u{q2). 

Proof : First we observe that (3.16) imphes the coefficients of u satisfy (3.11). 
This is obvious for c^i and c^s, but it follows for Ck2 because dnf{qo) — J hfd/j, for 
a fixed harmonic function h. Now apply Lemma 4.2 to the function u — u to obtain 



u{qi) - u{qi)\ = J VkA^{u - u)diJ, 



< cR'^\\vk\ 



But we easily obtain ||vfe||oo = 0{X2^) from (4.6) and the conjectures. Letting 
A; — > oo we obtain u{qi) — u{qi) = 0. Q.E.D. 

Proof of Theorem 4.1 : We begin by proving u = u under the assumption that 
u is even and R < Af*. Since A^u satisfies the same hypotheses as u, we conclude 
from Lemma 4.4 that A^[u — u) vanishes at all three boundary points, for any j. 
Let G{x, y) denote the Green's function and G^{x,y) the j-fold iteration of G. The 
vanishing at boundary points means that 

(4.7) u{x) - u{x) = J &{x,y)A^{u{y) - u{y))dfi{y). 
We have an explicit representation 

oo 

(4.8) G^{x,y) = J2{^^r'M^)My) 

k=l 
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for an orthonormal basis of Dirichlet eigenfunctions {(pk} with —Aipj^ — X^ipk- 
This yields the estimate 

(4.9) ( [ I \GKx,y)\^d^{x)d^{y)Y' = < c{\?r^ 

•^•^ fc=i 

by the Weyl asymptotics of {A^}. Thus 

\\u-u\\2 < c{X^)~'\\A^{u-u)\\2 < c{\^)-^W. 

Letting j — > oo we obtain ||tt — •u||2 = hence u = u a& desired. 

Next we can remove the assumption that u be even by writing m as a sum of 
even functions about each of the three boundary points using (2.26). It is clear 
that the hypotheses on u are inherited by the three summands, and a sum of three 
entire analytic functions is entire analytic. 

Finally, we need to relax the assumption that R < Af to i? < A2. To do this 
we consider u o for all words of length 2 (because 5~^A2 < A^). Then u o Fy^ 
satisfies (3.16) with R < Af, so by the previous argument it is entire analytic. 
This means for each w there exists entire analytic with u = on Fu,{SG). 
Next we claim that uqo = ttoi = uo2- To see this we may assume without loss 
of generality that uqq = by replacing u hy u — uqq. So u is assumed to vanish 
on Fq{SG), and we need to show that it vanishes on Fo{SG). By Lemma 4.4 we 
have u{Foqi) — u{Foq2) — 0, and more generally A^u{Foqi) — A^u{Foq2) = by 
the same reasoning for A^u. Let us consider uqi which equals u on FoFi{SG). 
At the point F^qi where the cells FqFi{SG) and Fq{SG) intersect, we have A^u 
vanishing and also d^A^u vanishing (obvious for the normal derivative with respect 
to Fq{SG), and then true with respect to FoFi{SG) by the matching condition for 
normal derivatives). Thus the local power series expansion in FqFi{SG) of uqi 
about the point F^qi contains only Pjs terms, so uqi and more generally A^uqi 
must be odd, so the vanishing of A-'ttoi at the second boundary point Foqi implies 
the vanishing at the third boundary point FoFiq2. So our previous argument shows 
that Uqi is identically zero. 

The same argument works in the other two cells of level one, so we now know 
that there exist entire analytic functions uq, ui, U2 such that u = Uj on Fj{SG). 
We need to show uq = ui = U2, and by subtracting uq we may assume without loss 
of generality that uq = 0. At this point we cannot simply repeat the argument of 
the previous paragraph because the cell Fi{SG) is too big. Of course we can argue 
as before that Ui and more generally A^ui vanishes on all three boundary points 
of Fi{SG), and that it is odd about the vertex Fq^i. It is this oddness that saves 
the argument. Instead of (4.7) for ui o Fi we have 



(4.10) 



«ioFi(x)= / G^{x,y)A\uioFi){y)dfi{y) 
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where & denotes the j-fold iteration of the odd part of the Green's function. 
Instead of (4.8), & has the same representation where the sum is restricted to the 
odd eigenfunctions. The eigenfunction associated to Af is even, so the smallest 
eigenvalue appearing is X2 ~ 55.8858 Thus we obtain the estimate 

\\ui o Fi\\2 < c{X^)-^5-^ W , 

and this shows ui = because A2 < 5A^. Q.E.D. 

It is interesting that the growth conditions (3.16) imply the specific identities 
(2.22). There is nothing analogous to this in the theory of real analytic functions. 
In some way it is reminiscent of the Cauchy integral formula for complex analytic 
functions. But we don't want to read too much into this, since (2.22) holds for 
nonanalytic functions as well. 

Corollary 4.5. If u is defined on a cell Fyj{SG) and satisfies 
(4.11) ||A%||ioo(p^(5G)) = 0{R^) for some R < X2 

then u has a unique extension to an entire analytic function. 

Proof : The theorem shows uoFyj is entire analytic. Then apply Theorem 3.6. Q.E.D. 

We can also consider entire analytic functions on any infinite blow-up of SG. 
The coefficients must satisfy (3.11) for all > 0, and the characterization requires 
the estimate (3.16) to hold locally for all R> 0. 

§5. Expansions about junction points. 

A junction point is a boundary point of two cells, so an entire analytic function 
will have two different local power series (3.14) centered at the point, each valid 
in a different cell. Since each local power series determines the function, it also 
determines the other local power series. Since the coefficients of the local power 
series are just the jets at the point with respect to each cell, these jets determine 
each other. The first goal of this seciton is to make this determination explicit. 

To be specific, consider the junction point Fogi = Fiqo. We will write Fqqi = qoi 
and write 

(5.1) {A^u{qoi),dnA^u{qoi),dTA^u{qoi)) 

for the jet associated with the cell Fo{SG), and FiQq = qiq and 

(5.2) {A^u{qio),dnA^u{qio),dTA^u{q^o)) 

for the jet associated with the cell Fi{SG). We know some relationships between 
the jets (5.1) and (5.2), namely 

(5.3) A^u{qoi) = A-'it(gio) and dnA^u{qoi) = -dnA^u{qio). 
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Note that (5.3) is valid for all u e domA°°, but there should be no connections 
between tangential derivatives without the assumption that u is an entire analytic 
function. On the other hand, for entire analytic functions, we expect an identity of 
the form 

oo 

(5.4) dTu{qoi) + dTu{qio) = ^ YidnA^u{qoi) 

to hold for certain coefficients Y^. Note that (5.4) applied to A^u yields 

oo 

(5.5) dTA^u{qoi)+dTA^uiqio) =Y,^i-j9n^''^iqoi), 

e=j 

and (5.3) and (5.5) show how the jets (5.1) and (5.2) determine each other. We 
may also interpret (5.4) as a matching condition for tangential derivatives. 

Our strategy for determining the Y coefficients will be to first consider the case 
when M is a polynomial, making the sum finite. It is convenient to consider the 

(2) 

monomials Pjj^ , because the p2 symmetry is also a symmetry about qoi. For even 
functions, both sides of (5.4) are zero regardless of the Y coefficients: the left side 
vanishes because of the oddness of the tangential derivative, and the right side 
because of the matching condition dnA^u{qoi) = —dnA^u{qoi) and the evenness 
of the normal derivative and Laplacian. Thus we need only check (5.4) for the 
monomials -Pjl"*- 

Lemma 5.1. The rriMching condition (5.4) holds for all polynomials for the Y 
coefficients satisfying Yq = 4 and recursively 



(5.6) Yj = -aj - 18j2^j+i-£-^ + Zl^^((2 ~ 

e=o £=0 

+ ^(5Q!j+i_^_fe)nfe5 - 3nj+i-i-kak5 for j > 1. 



fe=o 



(2) (2) 

Proof : When j = we compute directly that OtPqs \Qqi) + drPo^ (qio) = —4 

(2) 

and dnPgs (qoi) = —1, so Yq = 4. For j > 1 we use Corollary 2.14 to rearrange 
Pj^ around go- By (2.43) we obtain 

(5.7) P^!^ = -^P}? + 3 ^(a,+i_,Pi? + n,+i_,4°)). 



(2) 

Because P^^^ is odd we have 



dTPS\qoi) + dTP\l\qio) = 2dTP\l\qoi). 
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By (5.7) and Theorem 2.12 we have 

3 

(5.8) 2dTP^l\qoi) = aj + 18 ^ 



'5^ 



and 



(5.9) dnPjl\qoi) = ~ + X](5aj+i-fe'^fe5 - 3nj+i_feafe5 ^). 



fe=o 



amce i^^^ - 



(2) 

we have that (5.4) for u — P^^ yields 



Y, = 5]F,a„p(;!,)3(goi) - 2dTPg\qoi). 

£=0 

Substituting (5.8) and (5.9) yields (5.6). 

Conjecture 5.2. The coefficients Yj satisfy 
(5.10) \Yj\<cX2^. 

The numerical evidence for Conjecture 5.2 is presented in Table 5.1. 



Q.E.D. 



J 

1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



-4. 

-0.08888888889 
0.0002304526749 

-0.1434871749 10"^ 

0.1023938272 10"^ 
-0.7503519662 lO'^^ 

0.5527533783 lO-^^ 
-0.4076138308 lO'^^ 

0.3006465014 lO'^^ 
-0.2217590148 lO'^^ 

0.1635723837 lO'^o 
-0.1206533528 lO-^^ 

0.8899568485 10-^^ 
-0.6564452839 lO'^^ 

0.4842037197 lO'^^ 
-0.3571558034 lO'^^ 

0.2634433871 lO-^^ 
-0.1943197270 10"^^ 

0.1433330961 lO'^^ 
-0.1057246052 lO'^^ 

0.7798402782 lO'^^ 

Table 5.1 



{-x2yYj 

-4. 
12.05085573 
4.235674447 

3.575397353 
3.459038654 
3.436505741 
3.432052039 
3.431166398 
3.430989845 
3.430954602 
3.430947563 
3.430946155 
3.430945874 
3.430945818 
3.430945807 
3.430945805 
3.430945805 
3.430945805 
3.430945805 
3.430945805 
3.430945805 
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Theorem 5.3. Assume Conjecture 5.2. If u is any entire analytic function, then 
(5.4) and (5.5) hold for the Y coefficients given in Lemma 5.1. More generally, if 
X is any junction point in Vm+i \ Vm, then 

00 

(5.11) drA^uix) + d^A'u{x) = J] 3'^5-"'^^-'^Ye-jdnA^u{x), 

t=j 

where Ot and are derivatives with respect to the left cell at x and is the 
derivative with respect to the right cell. 

Proof : Note that the right side of (5.4) converges absolutely. The issue is then 
whether the tcrni-by-term differentiation of power series extends to normal and 
tangential derivatives at points other than the expansion point. For normal deriva- 
tives this is easy to see because of the integral representation. But in any case this 
follows by combining Theorem 3.4 (the explicit expression (3.12) for the rearranged 
coefficients) with Theorem 3.1 (the jet formula (3.5) at the expansion point). We 
then obtain (3.10) by applying (3.5) to the function no for \w\ = m. Q.E.D. 

Next we consider the question of what would be a natural notion of a power series 
expansion centered about a junction point. We will see that there is no completely 
satisfactory answer. Again to be specific we consider the point qoi = qiq. We would 
like to have at least the following four conditions holding: 

(i) every entire analytic function has an expansion; 

(ii) the expansion is valid in a neighborhood of goi, perhaps Fo{SG) U Fi{SG); 

(iii) the individual terms are polynomials that vanish to higher and higher order 

near goi; 

(iv) the rate of growth of the coefficients should be characterized for entire ana- 
lytic functions. 

The local power series with respect to one of the cells, say Fo{SG), gives a 
satisfactory answer only on that cell, but if we continue those monomials around we 
will find that the vanishing rate near gio is not satisfactory. In fact the tangential 
derivatives will have to be nonzero by Lemma 5.1. For this reason we consider 
carefully what it takes to meet condition (iii). We denote by Pjl^^ the monomials 
of the Fq{SG) local power series about qqi, so that 

dnA'P^l'\qoi) = SjeSk2 

dTA'P^l'\qoi)^d,idk3 

or more precisely 

P^','\x)^5-^P^I\f-'x) 
<)(.) = ^5-P,«(Fo-M 
P^','\x)^5-^-'pff{F,-'x). 
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Note that Pj°^^ and Pjg^^ extend to even polynomials about goi: so they will have 

the same vanishing rate on both cells. We want to replace -Pj2^'* by a different 

polynomial Pj2'^ that will have the same j-jet (except for dx u{qoi)) ^ but will 
extend to be odd. This will give it the correct order of vanishing, but in exchange 
we have to take a higher order polynomial. The lowest possible order is 2j: 

^ j 

£=0 

for the appropriate choice of constants. Note that we can exclude -P^^^^^^^ terms be- 
cause we want the possibility of odd extension. We will take ajj = 1 in order to ob- 
tain the correct j-jet. The odd extension means dTA'^Pj2^^ (qoi) = 9t A"Pj2^'' (gio), 
so we have 2j + 1 equations of the form (5.5) to satisfy, and these will determine 
the remaining 2j + 1 constants. The equations are 

2j 

(5.13) 2aTA-pjf )(gio) = J] n_na,A^pj5^)(go2), 

fe=n 

and when < n < j the left side is zero and we obtain 

= X;^fe-n5nA'=Pj;'^(g0l) = J]Ffe-na,(2,-fe) 



k=n k—j 



SO 



J 

(5.14) = J2y^j-e-naj£. 

£=0 

We use these equations to solve for aj£. When n < j < 2j the left side of (5.13) is 
2&j(2j-n) so 

k=n 

and by letting £ = 2j — n we have 

1 ^ 

(5.15) bji = 2 X] ^kaj{£-k) for < £ < j. 

k=0 

In Table 5.2 we show the values of aji and bj£ for small values of j. It is difficult 
to discern a pattern in these results. We have obtained graphs of Pj2^'' for small 
values of j using (5.12), but it appears that round-off error becomes significant 
before any pattern emerges, so we are not able to offer any conjectures about the 
growth rate of these functions as j — > oo. 
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j I ajl bji 

1. 2. 

1 0.02252966406 0.04505932812 

1 1 1. 1.999249011 

2 6461.417615 12922.83523 
2 1 -39.86777272 -295.1161326 

2 2 1. 9563.195714 

3 0.1631072895 lO'^ 0.3262145790 lO'^ 
3 1 48581.69671 42794.29693 

3 2 -109.6002902 0.2411384099 lO'^ 

3 3 1. 86782.07999 

4 -0.1623039023 10^° -0.3246078045 10^° 
4 1 -0.6442287860 10^ -0.7474445645 10^ 
4 2 -299734.8354 -0.2399788368 10^° 
4 3 -347.4611669 -0.1101312661 10^ 

4 4 1. -751724.7199 

5 0.1010368178 lO^^ 0.2020736356 lO^^ 
5 1 0.4380632964 10^^ 0.5393372002 10^^ 
5 2 0.3374174349 10^° 0.1494085527 10^^ 
5 3 0.1015644445 10^ 0.7403235769 lO^^ 
5 4 -909.3198857 0.7249040413 lO^o 

5 5 1. 0.1921254540 10^ 

6 -0.1389829261 10^^ -0.2779658521 10^^ 
6 1 -0.6247328496 W^^' -0.7861892790 10^^ 
6 2 -0.5605362673 lO^^ -0.2055333917 lO^^ 
6 3 -0.2151475440 lO^^ -0.1051115464 lO^'^ 
6 4 -0.2169919676 10^ -0.1159983908 lO^^ 
6 5 -1787.130925 -0.4257054009 lO^^ 
6 6 1. -0.4383706038 10^ 
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i 


7 

t 




bji 


7 





0.1330959781 10^ 


0.2661919562 10 


7 


1 


f\ c\ A'\ c\'\ nc\cf\ ir\21 


n TO /I Tone o 1 ^7 1 r\21 

0.7847295317 10 


7 


2 


U.oU547aDo57 lU^ 


n 1 n£?oo£?CTi o t r\23 

0.19o83d5718 10 


7 


o 

3 


U.z7U75Uo9o7 lU 


0.1030137030 10 


7 


4 


0.4581523610 10 


0.1231428577 10 


7 


5 


0.2620127789 10 


0.5414059059 10 


7 


6 


3880.162356 


0.yi582oo227 10 


7 


7 


1. 


0.5243561927 10 


o 

8 





— 0.28493d7d»o 10 


A rr ^3 AOTO CT OTCT 1 a25 

— (J.5d987o5o75 1(J 


8 


1 


—0.1352804490 10 


0.1755939762 10 


8 


2 


—0.1478090302 10 


A /101/11T/1ACA 1a25 

— U.4z14174(jd9 1(J 


8 


3 


—0.7540725789 10 


—0.2261525660 10 


o 

8 


4 


— U.17o(Jod15od lU 


n ononci£;nT£; 1 n22 

—0.2930516976 10 


8 


5 


—0.1895987908 10 


A ir'iioiAr'iA iA2n 

—0.1511819510 10 


o 

8 


D 


— 0.775DD75150 10 " 


— 0.3520528934 10 ' 


o 

8 


7 


—3618.462380 


A O^AATOAOTA 1 a1 4 

— 0.3790729379 10 


8 


8 


1. 


—0.1552676258 10 


9 





0.4817483229 10 


0.9634966458 10 


9 


1 


0.2289760048 10^8 


0.2973692352 10^8 


9 


2 


0.2513117964 10^6 


0.7125008828 10^9 


9 


3 


0.1299020030 10^4 


0.3827224251 10^8 


9 


4 


0.3154544064 lO^i 


0.4977745865 10^6 


9 


5 


0.3859718201 lO^^ 


0.2600348690 10^^ 


9 


6 


0.2325380299 lO^^ 


0.6310751388 lO^i 


9 


7 


0.5539946952 lO" 


0.7717084596 10^8 


9 


8 


-6592.977986 


0.4652032965 10^^ 


9 


9 


1. 


0.1107495241 10^^ 



Table 5.2 
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§6. Exponentials. 

Eigenfunctions of the Laplacian give us a natural class of special functions on 
SG. Until now, most attention has been paid to eigenfunctions satisfying Dirichlet 
or Neumann boundary conditions, which forces the eigenvalue to be positive. In 
contrast, we will mainly explore negative eigenvalues in this section, so we are ex- 
ploring the analog of the functions cosh v^At and sinh ^/Xt on the unit interval and 
their extension to the positive real line. Of particular interest is the linear combina- 
tion that yields e~^*, the unique choice that exhibits exponential decay (either as 
A — > oo or as t — > oo) as opposed to exponential growth. It is embarrassing to note 
that the exponential e^* does not distinguish itself among linear combinations of 
cosh \/Xt and sinh ^/Xt, if one is forbidden to use odd order derivatives. So we have 
not been able to find its analog on SG. 

The space of all eigenfunctions with a fixed eigenvalue has dimension three, as 
long as one avoids Dirichlet eigenvalues. For fixed A > we can choose a basis C\, 
Sx, Q\ for the space of solutions to 

(6.1) -A« = -Xu 

determined by the conditions that C\ and S\ are even and Q\ is odd with respect 
to po) and 

(6.2) CM = 1, dnCxiqo) = 



(6.3) Sx{qo) = 0, dnSxiqo) = 

(6.4) drQxiqo) = 1 

where the normalization factor ax will be chosen later. This means that we have 
global power series representation 

oo 

(6.5) Cx{x) = J2X^P^?{x) 
and 

oo 

(6.6) Qx{x) = J2X^P^^\x), 

i=o 

and a local power series representation 



(6.7) 



oo 

Sx{x)^axJ2X^P^',\x) 

j=0 
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valid on Fq{SG) provided A < 5^X2- We may also use (6.5) and (6.6) on the 
blowups Fq'^{SG) for any n. Of course, none of these functions are entire analytic 
for A > A2. 

00 

We will consider the infinite blowup SGoo = [J Fq'^{SG) to play the role 

n=0 

of the positive reals vis-a-vis the unit interval. Of course there are uncountably 
many infinite blow-ups of SG. We have chosen the simplest one to study first. To 
understand the "behavior at infinity" of these functions it suffices to study the 
values at the points Xn = FqQi as n ^ —00, for we may then get the values at the 
points i/n = FQq2 by parity, and then fill in by spectral decimation. 

For S'Goo we have graphs for any integer n. Since —A is negative we never 
encounter the exceptional eigenvalues 2, 5 and 6. Thus the method of spectral 
decimation says that u satisfies (6.1) on SG^o if and only if the restriction of u to 
r„ is a graph eigenfunction with eigenvalue An, where {An}nez is a sequence of 
negative numbers characterized by 

(6.8) A„_i = A„(5 - A„) 
and 

(6.9) -A= lim ^5"A„. 

n— >oo 2 

Note that A^ — > as n — > 00 and An — > —00 as n — > —00. It is easy to see that the 
sequence {Aj} is uniquely characterized by these conditions, and the values may be 
effectively computed to any desired accuracy by replacing the limit in (6.9) by the 
value for a fixed large n and then using (6.8) to run n down. 

The fact that u restricted to F^ is a A„-eigenfunction means that if we take any 
cell of level n — 1 with boundary points a, 6, c, and if d is the midpoint between a 
and 6, then 

(6.10) u(i) = (^-An)(«(a) + «W)+2u(c) 

(2 - A„)(5 - Xn) 

(see [DSV] Algorithm 2.4). 

Lemma 6.1. The recurrence relations 

(a.n) ,,.,.(l-^±£zi^ 

and 

(6.13) Qx{xn) = 

O — An 

hold for all integers n. 

Proof : Apply (6.10) for a = go, b = F^-\qi), c = F^-\q2) andd = F^{qi). Q.E.D. 
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Lemma 6.2. The function Cx is positive. The function Sx, with the appropriate 
choice of ax, is positive everywhere except at qo where it vanishes. The function 
Qx vanishes on the symmetry line through qo and is positive on the qi half of the 
symmetry line. 

Proof : Because A„ < for all n, the coefficients in (6.10-6.13) are all positive. 
That means that if u is nonnegative on the boundary of a cell and strictly positive 
at one of the boundary points then it is strictly positive in the interior. Thus it 
suffices to show that Cx{xn)^ Sx{xn) and Qx{xn) are positive. For Sx and Qx it 
suffices to show Sx{qi) and Qx{qi) are positive, since we can solve (6.12) and (6.13) 
for Sx{xn-i) and Qx{xn-i) with positive coefficients. But we can make Sx{qi) > 
by the appropriate choice of sign (negative) for ax, and Qx{qi) > follows easily 
from drQxiQo) — 1- When we solve (6.11) we obtain 

r ^ {2-Xn){5-Xn)Cx{Xn)-{^-Xn) 

(6.14) Cx{Xn-i) = 7 , 

D — An 

which contains a negative coefficient. Nevertheless, if Cx{xn) > 1 then (6.14) 
implies 

r u W (2-An)(5-An)-(4-An) 

Cx{Xn-l) > > 1, 

6 - An 

SO it suffices to show Cx{qi) > 1. This follows because the contrary assumption 
Cxiqi) < 1 and (6.13) would imply dnCx{qo) > 0. Q.E.D. 

Theorem 6.3. (a) For all n we have 

(6.15) Cx{xn) = l-^. 
(b ) For the appropriate choice of ax we have 

Xn 



(6.16) Sx{xr,) = -^ll{l + 



fe=0 2-An-fe 

and hence 

(6.17) hm Sx{Xn)/Cx{Xn) = l. 

n— >— oo 

(c) For all n < we have 

(6.18) Q^(^^) = _3^ 
and hence 

(6.19) lim QxiXn)/CxiXn) = ^. 

n— > — oo A 
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Proof : (a) A direct calculation using (6.8) shows that 1 — ^ satisfies the same 
recurrence relation (6.11) as C\{xn). Thus if we define Cx{xn) ~ 1 — C\{c[o) = 1 
and extend C\ to all of SGoo using (6.10), we will have an even A-eigenfunction. 
But a direct computation shows 

■5yl 

3- 



dnCxiqo)^ lim (D'Ia, =0 



because \j = 0{5~^) as j oo. So C\ = C\, proving (6.15). 

(b) First we observe that the infinite product in (6.16) converges, because of the 
rapid growth of A^j, as n ^ — oo. Since (6.12) may be written (using (6.8)) 

(6.20) M^ = (i+ 4 ^5.(x_l) 



An ^ 2 — An ^ An-1 

it follows that the right side of (6.16) satisfies (6.12). Since S\ was only defined up 
to a multiplicative constant, we may choose a\ to make (6.16) hold. Note that from 

(6.20) we obtain S\{xn) = 0((|) ) as n — >• oo, which is consistent with Sx{qo) = 
and dnSx{qo) ^ 0. Then (6.17) follows from (6.15) and (6.16) by inspection. 

(c) We may rewrite (6.13) as 

QxjXn) _ QxjXn-l) 
An An-1 

using (6.8), hence Qx{xn) = XnQ\{xo) for all n. But then 

1 = drQxiqo) = Hm 5"(Q(xn) - QxiVn)) 

n — >cxD 

= 2Qx{xq) lim 5" An 

n— >oo 

4 

This proves (6.18), and then (6.19) follows by inspection. Q.E.D. 

Wc can compute the value of ax = dnSx{qo) exactly. From the definition and 
(6.16) we have 

/5 \ 

(6.21) dnSxiqo) = -'2 lim - Sx{xn) 

n— »oo \6 / 



lim ^(-Y TT (l ^ 

fc=0 ^ 



= -1a lim -ifTri+ t ) 

3 n^<^3-llV 2-An-fc/ 
fe=0 

=-i^n(i-jz^)„'-n(|f^) 
=-iAn(^+^)n(S). 



j=0 " ^' k=l 
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Definition 6.4 : For A < define the decaying exponential function E\ by 

(6.22) Ex{x) = Cx{x)-Sx{x). 
Theorem 6.5. Ex{xn) = 0{X-^) as n ^ — oo. In fact 

(6.23) Urn Xr^Ex{xn) = -l 

n—^—oo 

and 

(6.24) Urn Cx{xnf-Sx{xnf = ^. 
More precisely 

(6.25) Exixr.) = -\- + + , X ^ + ^(^n')- 

Z — An Z — An-l [Z — An)[^ — An-l) 

Proof : From (6.16) we obtain 

(6.26) 5.(.„) = (l + ^) (l + + 0(A„-3) 

because Xn/Xn-2 = 0{X~^). Substituting (6.26) into (6.22) and using (6.15) we 
obtain (6.25). Using (6.8) we see that the first two terms on the right side of (6.25) 
sum to 



2 — Xn 2 — 5A^ + Xj 

The third term is clearly 0(A~^), so we obtain (6.23). From (6.26) we find Sx{xn) = 
0(1) and this yields (6.24). Q.E.D. 

Note that (6.26) and (6.25) allow for the efficient computation of Sx and Ex for 
n sufficiently negative. On the other hand (6.22) is computationally unstable since 
it involves subtracting values that are large and nearly identical. In Table 6.1 we 
present some numerical computations of these functions. 

Instead of fixing A and taking the limit as n — > — oo, we could look at values 
at xq and let A — > — oo. As long as |Ao| is large, (6.25) and (6.26) will be good 
estimates. Table 6.2 shows this behavior. We could also allow A to be complex, as 
long as the real part is positive to avoid the exceptional values for An- 

We now turn our attention to eigenfunctions with positive eigenvalues, with the 
goal of using information gleaned from spectral decimation to shed some light on 
the recursion relations from Section 2. Keeping the same notation as before, we are 
interested in the function 

oo 



C-xix)^Y^{-xypj,{x) 

i=o 
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-.2923335134 10^^ 


.7308337835 10^^ 


.7308337835 10^^ 


-5 


-.8545888306 10^^ 


.2136472076 lO^^ 


.2136472076 lO^^ 


—6 


— .7303220694 10'^ 


.1825805173 10'" 


.1825805173 10'" 
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Q 

—6 


.O/oyZdDOOO lU 


.io4yoz/Uoy iU 


— .yyyyyyyy4o 


A 

—4 


.2048759594 10 


.3420750458 10 


— 1.0000000000 


-5 


.5989210902 lO^^^ 


.1170153370 10-34 


-1.0000000000 


-6 


.5118312741 10^^ 


.1369258909 10"^^ 


-1.0000000000 


-7 


.3738016753 10^^^ 


.1874869960 lO-^^g 


-1.0000000000 


-8 


.1993747210 lO^^s 


.3515137367 10'^^^ 


-1.0000000000 


-9 


.5671889891 10^^^ 


.1235619071 10-^^s 


-1.0000000000 


-10 


.4590322393 10^^^^ 


.1526754489 lO-^^^^ 


-1.0000000000 



Table 6.1. Values of functions at X-j for A = 10.70160380. 



and its values at the special points xq — qi and xi = Foqi. It is convenient to define 
A„ (here we only care about n > 0) to satisfy (6.8) but to remove the minus sign in 
(6.9). For the Dirichlet and Neumann eigenfunctions we know exactly what these 
values are, and then we can use Theorem 6.3 (a) to conclude that C-\{xo) = 1 — ^ 
and C-x{xi) = 1 — (Strictly speaking, we need to use an analytic continuation 
and limit argument to get this for the values we are interested in.) In particular, if 
Ao = —6 then C-x{xo) — 5/2, or 



oo oo 

5^(-A)^P,.i(gi) = 5](-A)^a,=5/2. 

i=o j=o 



CALCULUS ON THE SIERPINSKI GASKET 



49 



Ao 



A 



Ex{xo) 



-100 44.19536761 

-500 87.71437197 

-1000 112.0105482 

-5000 182.0354932 

-10000 218.2833208 

-50000 317.2473555 



.009711493217 
.001988095160 
.0009970119472 
.0001998800959 
.00009997001199 
.00001999880010 



first 2 terms 
in (6.25) 

.01008584733 
.002003881410 
.001000985089 
.0002000398801 
.0001000099850 
.00002000039988 



first 3 terms 
in (6.25) 

.009712435727 
.001988103065 
.0009970129413 
.0001998801039 
.00009997001299 
.00001999880010 



Ao 



A 



Sx{xo) 



first 2 factors 
in (6.26) 



first 3 factors 
in (6.26) 



-100 


44.19536761 


25.99028851 


25.98039216 


25.99028756 


-500 


87.71437197 


125.9980119 


125.9960159 


125.9980119 


-1000 


112.0105482 


250.9990030 


250.9980040 


250.9990030 


-5000 


182.0354932 


1250.999800 


1250.999600 


1250.999800 


-10000 


218.2833208 


2500.999900 


2500.999800 


2500.999900 


-50000 


317.2473555 


12500.99998 


12500.99996 


12500.99998 



Table 6.2. Values of functions at xq for various A values. 

This happens when A = A2, the second nonzero Neumann eigenvalue (not to be 
confused with the A2 in (6.8) and (6.9)). This allows us to compute the limit of 
Pj/tj+i as j — > 00. Indeed, from (2.34) we have 



A. 
tj+i 

We expect to have 



= 6 ^ aj+i. 

£=0 



t£ 

£=0 



'^j+l-£ 
tj+i 



6. 



tj + l-£ 
tj+l 



and so 



lim 



6 ai{-X2Y - 6 = 6- ^- 6 = 9. 



This is confirmed by the data in Table 2.2. 

We are also interested in the solutions of the equation 



(6.27) E«H-^)' = -2- 

£=0 

This holds for z = A2/5, because in this case Ai = 6, and 

00 

C_a(xi) = J]a,(-A2/5)^. 

£=0 
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But it also holds for 2; = Af, because in this case Aq = 6. In fact it is easy to see 
that Af is the smallest solution of (6.27) (there are infinitely many other choices of 
A with either Ai = 6 or Aq = 6). Figure 6.1 shows the values on Vi of the function 
C_A in these cases. 



Figure 6.1 : The values of C-x{x) on Vi vertices for (a) Aq = —6 and Ai = 6, 
(b) Aq = 6 and Ai = 2, (c) Aq = 6 and Ai = 3. 

We can now explain why the recursion relation (2.11) for (3j is unstable. It is 
clear by inspection that the middle term on the right side of (2.11) is much larger 
than the other terms, so we would expect that a solution of (2.11) would be close 
to a solution of 

2 

/5. = ^:tE«^-^5^"'/^^, 



£=0 



which may be rewritten as 



(6.28) -^ = E«^5-^%^. 

If we look for a solution of (6.28) of the form Pj = {—5z)~^ then we obtain 
^ 1 

'^^a£{—zY — — , which is very close to (6.27) in view of the very rapid de- 

cay of ai. The solution to (6.28) should thus be an infinite linear combination of 
exponential solutions with z a solution to (6.27). In the generic case the dominant 
term should correspond to the smallest solution of (6.27). Thus we expect the solu- 
tion to (6.28) to behave like a multiple of {—5Xf)"\ and numerical computations 
confirm this. This pseudo-solution of (2.11) attracts any approximate solution of 

(2.11) that strays from the exact solution. 

00 

A related observation is that ^ a£{—zY = 1 holds for 2; = Af ~ 55.885828 . . . 

e=o 

00 

by (6.15), since in this case Aq = and Ai = 5. In the form ^a£(— Af )^ = this 

i=i 

suggests that the entries of the matrix a{a)~^, which are just 6Tj, should decay 
like (— A|^)~-'. The numerical data in Table 6.3 confirms this. This explains the 
instability in the recursion relation for {tj}. 
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j 









1. 


1. 


1 


- 03333333333 


1 862860915 

_L « \^ vy Kj vy vy t_y -l 


2 


0007407407407 


2 313500526 

« _i_ vy vy vy vy ^ vy 


3 


- 00001433691756 


2.502423700 


4 


2637965601 10~^ 


2 573213790 


5 


- 4766054541 10"^ 


2 598169232 

• \j %_j v_y -i- vy \j ?_y 


6 


8556101104 10~^° 


2 606669803 

^< vyvyvyvyvyty w vy vy 


7 


- 1532663873 10"^^ 


2 609508520 

^ • vy vy ty vy vy vy vy ^ vy 


8 


2743475872 10"^^ 


2.610445492 


9 


- 4909650195 10~^^ 


2 610752605 


10 


8785480907 10"^'^ 


2 610852844 


11 


- 1572060595 10"-^^ 

• j-ty 1 ^vyvyvy^yty^y j. vy 


2 610885478 

^ • vy -1- vy vy ^ i vy 


12 


2812997595 10~^° 

• i-j w -i-^tytyi tyt_yty j. vy 


2 610896085 

^ • vy -L vy vy ty vy vy vy vy 


13 


-.5033478852 lO'^^ 


2.610899530 


14 


.9006721805 lO-^^ 


2.610900647 


15 


-.1611629185 10-25 


2.610901010 


16 


.2883788845 10"^^ 


2.610901127 


17 


-.5160143489 lO'^^ 


2.610901165 


18 


.9233366935 lO'^i 


2.610901177 


19 


-.1652183992 lO-^^ 


2.610901182 


20 


.2956355963 lO-^^ 


2.610901182 




Table 6.3 





We also observe that the values of C-x^{x) given in Figure 6.1 (a) show that the 
rearranged power series at qi does not converge to C-x^ outside the cell Fi{SG). 
Indeed, the even part of the power series about qi, if it converged in SG, would have 
to be f Ei->^2yP^l\x), which gives the incorrect value of 25/4 for |(C_a2(9o) + 
C-A.(g2)) = 7/4. 
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